Dynamically localized systems: entanglement exponential sensitivity and efficient 

quantum simulations. 
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We study the pairwise entanglement present in a quantum computer that simulates a dynam- 
ically localized system. We show that the concurrence is exponentially sensitive to changes in 
the Hamiltonian of the simulated system. Moreover, concurrence is exponentially sensitive to the 
"logic" position of the qubits chosen. These sensitivities could be experimentally checked efficiently 
by means of quantum simulations with less than ten qubits. We also show that the feasibility of 
efficient quantum simulations is deeply connected to the dynamical regime of the simulated system. 

PACS numbers: 05.45.Mt, 03.67.Mn, 03.67.Lx, 07.05.Tp 



The study of the entanglement in quantum chaotic sys- 
tems is object of growing interests These works are 
based on classical simulations of quantum systems. De- 
veloping quantum computers will enhance our investiga- 
tions power giving the possibility of performing nowadays 
unaccessible simulations. However, being able to perform 
exact quantum computation is not a sufficient condition 
to assure that quantum simulations will be useful: ex- 
tracting useful information in an efficient way is a dif- 
ficult task that, in some cases, drastically decreases the 
efficiency of quantum algorithms @,0|. Moreover, in Q| 
it has been shown that the entanglement present in a 
quantum system plays an important role in the possibil- 
ity of performing efficient classical simulations of quan- 
tum systems. Thus, the presence of entanglement and 
the processes of quantum measurement may determine 
the efficiency of quantum algorithms. 
It is well known that quantum chaotic systems display ex- 
ponential sensitivity to small changes of the Hamiltonian 
[U in dynamical quantities as the fidelity [g ■ In it has 
been shown that the entanglement evolution is influenced 
by the presence of quantum chaos in a system following 
the characteristic decays of the fidelity. This result has 
been found in a system with no classical analogue. Here, 
we study the concurrence evolution between the qubits 
in a quantum computer that is running a quantum al- 
gorithm to simulate a quantum chaotic system, in the 
regime where dynamical localization occurs. Dynamical 
localization is another signature of quantum chaos: it has 
been first predicted in [lfj, and experimentally observed 
in cold atoms . It has been related to Anderson local- 
ization ^3 °f electrons in disordered crystals. Then, it is 
interesting to perform simulation on many-body or clas- 
sically chaotic quantum system to check whenever they 
localize or not: The analysis of localization length as a 
function of the parameters of the system is important to 
study transport phenomena. In |3( it has been shown 
that localization length can be extracted with polyno- 



mial gain if an efficient quantum algorithm to simulate 
the system exists. This is not a trivial property, as the 
efficient extracion of many interesting quantities is pre- 
cluded due to the exponential number of measurements 
needed to know all the wave function coefficients. 
In this paper we show that the presence of an exponen- 
tially localized wave function implies exponentially sensi- 
tivity in the concurrence between qubits. This sensitivity 
is both with respect to small changes in the Hamiltonian 
and to the pair of qubits chosen. The pairwise entan- 
glement can be extracted efficiently from any quantum 
algorithm that simulates a quantum system efficiently. 
Indeed, the measure of the reduced density matrix of 
two qubits is an efficient process that does not scale with 
the size of the system. Thus, both localization length 
and concurrence can be measured efficiently by means of 
quantum simulations. We also study the entanglement 
of blocks of qubits with respect of the other and we show 
that these properties are strictly related to the localiza- 
tion length. Thus, following |4(, the localization length 
determine the feasibility of efficient classical simulation 
of the system. 

The Quantum Sawtooth Map (QSM) is a suitable 
model of quantum chaotic system that displays dynam- 
ical localization which can be efficently simulated with 
a quantum computer. The QSM belongs to the kicked 
map family and it is defined by the Hamiltonian 



J25(t-rT), 
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where <5t(£) is a delta of Dirac and r £ Re, h and 9 
are canonical conjugate variables such that [6, h] = iT 
(we set h = 1, thus T = K e ft). Correspondingly, the 
quantum dynamics is described by the Floquet operator 
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where < 6 < 2tt and -N/2 < n < N/2. The wave 
function time evolution is computed by repeated appli- 
cations of operator (0) to the initial wave function, i.e. 
\ip(m)) = U m \il)(0)). The dynamics of the QSM is gov- 
erned by two parameters k, T, while the dynamics of the 
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FIG. 1: Concurrence of qubits 1 and 3 as a function of time for 
different map parameters K = y/2, k = K/T, T = 2-KM/2 71 " , 
n q = 10. From bottom to top M = 10 4 ,5 ■ 10 3 ,3 ■ 10 2 ,2 • 
10 3 , 5 • 10 2 , 10 3 , 7 • 10 2 and I = 0.03, 0.1, 31, 0.7, 11, 2.8, 5.7 re- 
spectively. Inset: peak of an exponentially localized wave 
function. The grey thick line follows the law fj. 

classical correspondent map depends only on the classi- 
cal parameter K = kT and displays chaotic dynamics 
for K > and K < -4 8]. For -4 < K < the map 
is described by mixed phase space (both integrable and 
chaotic). The classical limit of map is recovered for 
k — > oo and T — » keeping K constant. For T < 1 and 
K > 1 unique quantum feature appear and QSM displays 
dynamical localization. The classical chaotic diffusion in 
momentum is suppressed by quantum interference and 
the wave function is described by an exponentially local- 
ized wavefunction of the form 

\i>(n)\ « e-W l /Vi. (3) 

The inset of Fig^shows such typical localized wave func- 
tion. Any further evolution of the wave function is sup- 
pressed excepted from quantum fluctuations. The QSM 
in the localized regime is described by the time indepen- 
dent wave function @ with £ 7r 2 fc 2 /3. Dynamical lo- 
calization occurs after time t* ~ £ 9] . We stress the fact 
that the expression is a typical signature of quantum 
chaos. Notice that dynamical localization can be seen 
only if £ << N: In the thermodynamic limit dynamical 
localization always occurs. 

In |13| an efficient quantum algorithm has been pre- 
sented to simulate the map (|2J|. The algorithm is based 
on the quantum Fourier transform and it exploits the 
exponential efficiency of the quantum Fourier transform 
with respect to the classical fast Fourier algorithm [l4| . 
If performed on a rt g -qubits quantum computer, it needs 
0(riq) elementary quantum gates, while a classical com- 
puter needs 0(N — 2 nq ) elementary gates. Moreover, 
the QSM algorithm has no need of ancillary qubits: it is 
already possible to see dynamical localization with less 
than ten qubits. Notice that any quantum algorithm 
to simulate a quantum system is based on the introduc- 
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FIG. 2: Concurrence of qubits 1 and j as a function of time for 
different map parameters K — \/2, k = K/T, T = 2-kM /2 n " , 
n q = 10, M = 800, I = 4.36. From top to bottom j = 
2,3,4,5. 



tion of the coding of the dynamical variable in binary 
representation. In our case, the coding is defined by 
n = ') where n are the eigenvalues of h describing 

the QSM and ati = 0, 1 depending on which level of the 
i — th qubit is populated. Thus, the qubits are ordered as 
they represent different logical values, from the least sig- 
nificant qubit (i — 1) to the most significant one (i = n q ). 
Each of them naturally introduce a coarse graining of the 
space of the n |lfj| . From now on, we refer to the differ- 
ence of coarse graining level A(i,j) = 2 l — 2 J as coarse 
graining distance between the qubits. Moreover, for the 
sake of simplicity, we will identify the logic label i with 
the spatial position. However, our results are general as 
they can be recovered with an appropriate mapping be- 
tween spatial and logic position. It is common wisdom 
that entanglement is a necessary resource to exploit the 
quantum computational gain. It is then natural to study 
the pairwise entanglement in the qubits that describe a 
dynamically localized system. We show that pairwise 
entanglement is present between the qubits and that its 
value depends by the natural ordering introduced by the 
quantum algorithm. Indeed, differently from other stud- 
ies, the coding n = ^2i&i2 l , introduce a hyerarchy that 
was not present in spin systems as in 0. We show that 
in the regime of dynamical localization the concurrence 
value exponentially depends on the coding position of 
the pair of qubits and that it is exponentially sensitive of 
small change of the kick strength k. 

We quantify the pairwise entanglement present be- 
tween any pair i, j (with i < j ) of qubits by means of con- 
currence defined as C = max{Ai — A 2 — A 3 — A4, 0}, 
where the Aj's are the square roots of the eigenvalues 
of the matrix R = p 1 ^ , in descending order; the 
spin-flipped density matrix is defined by p 1 ' 3 = (a\ ® 
<rJ)(/9 lJ )*(of ® dj) (in this definition the two qubits ba- 
sis lajay) must be used). The reduced density matrix 
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FIG. 3: Concurrence saturation values between qubits i = 1 
and j as a function of localization length £. Different sim- 
bols represents different qubits couples: j — 2 (circles), 
j = 3 (squares), j = 4 (diamonds), j = 5 (triangles). 
Dashed lines follows the predicted behavior C(£) ~ £ and 
C(€) ~ exp(— j4€). Inset: Concurence values for £ = 2.8 and 
different qubits couples: i = 1 (circles) and j — i = I (squares). 
Full line represents C(j) ~ 2 J '. 



p*>J is given by 

n, l 

= II E ^a 1 ... A ...A»...a n ,Ax 1 .../9,.../S m ...a ni (4) 
k—l,k^i,j a k — 



0, 1 ... 3, the /3fc are the binary code of I and 

ft.../3 m ...Q„ = («1 • • • A • • • An • • I^W)- 



where i , m 
m and A Q1 

We start with an initial wave function equal to a mo- 
mentum eigenstate \ip(t = 0)) = |n), we compute the 
time evolution at a given time t by repeated application 
of the Floquet operator and at each time we compute 
the concurrence by means of the reduced density matrix 
Q. We plot in Fig^the concurrence evolution of two 
qubits as a function of time for different k values and 
fixed K: this choose of the parameter values corresponds 
to the same classical dynamics but different quantum dy- 
namics with different localization time and length. The 
stationary regime is reached for larger times as I value is 
increased. In Fig. |3we plot the concurrence evolution for 
fixed I and different qubit pairs. The different saturation 
levels and times are visible. Notice that the concurrence 
saturation level is not a monotonic function of the lo- 
calization length: indeed the two limiting situations, a 
momentum eigenstate and a random superposition of all 
the momentum eigenstate i.e. an ergodic wave function, 
have zero concurrence values for any qubit pair 0, thus 
a maximally entangled case should lie in between. 

In Fig. [3]we show in details the saturation value of the 
concurrence between the less significant bit (i = 1) and 
the others (j = 2,3.. .) as a function of the localization 
lcnght. There are two distinct regimes: for I < 1 the 
concurrence increases as C(£) ~ VI, while for I » 1 the 
concurrence is exponential sensitive to the localization 



length. The squared root dependence can be understood 
as follows: The number of addend different from zero 
in JU can be estimated as N e = £/Q = £/(A 2 j ). In- 
deed, the number of pairs of coefficients A ai _ p t _ _ _ an 
of coarse graining distance A(i,j) in an ensemble of £ co- 
efficients different from zero drops linearly with A(i,j), 
thus exponentially with j. The structure is repeated each 
2 J , giving the scaling Q ~ A V . For I < 1 only few levels 
of the wave function are significantly populated, and they 
are "near" , i.e. the binary code of their position n in the 

wave function (/3icv2 Pj ■ ■ ■ a n q ) differs only for the 

value of (3\. To estimate the sum we aproximate the 
wave function with the expression 



i), 



(5) 



where 4>i is a random phase and in + 1/2 the position 
of the peak of the wave function 0. Thus the density 
matrix will be composed by diagonal terms porportional 
to pW ~ 1/^. Notice that either p % Q q or p^ J 3 will be neg- 
ligible as they have a coarse graining distance is A >> I. 
Equivalently, the off diagonal matrix elements which dif- 
fer significantly from zero are those who are composed by 
products of terms which have the smaller coarse grain- 
ing distance, that is, the first lower and upper diagonal. 
Furthermore, if p§\ is negligible also p 1 ^ may be ap- 
proximate by zero. We then compute the matrix, 
which mix up the density matrix elements, leading to an 
approximate two by two matrix of the form: 



/0 0\ 

o ±# A o 1 
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The squared root of the non-zero eigenvalue of R p ma- 
trix is then proportional to Vi/G ~ Vi/2 2 ^ . Equation 
© implies that the saturation value of the concurrence 
for small I scales as C ~ Vi~/2 2: >: This prediction is in 
agreement with the numerical data presented in Fig. |3J 

For £ » 1 the concurrence has a completely different 
behavior (Fig. |3J): the number of coefficients different 
from zero in 10} increases, the approximation that leads 
to |[BJ| is not valid any more as the off-diagonal coeffi- 
cients of |@J| start to be not negligible. The overall effect 
is that the concurrence drops exponentially with £, up 
to a critical value where it drops exactly to zero. This 
happens when all the off diagonal term are on average 
equal. As £ oc k 2 , small changes in the Hamiltonian ||TJ 
of the kind k — » k+Ak change exponentially the pairwise 
entanglement present in the system. 

In the inset of Fig. [3]the concurrence saturation values 
for different qubit pairs are plotted for a given value of the 
localization length. As it can be clearly seen the points 
follows the estimate C(£) oc 2 _2j . In Fig. 0|we show that 
the two different regimes of the concurrence are typical of 
any pair of the qubits: For every couple of qubits chosen 



FIG. 4: Concurrence saturation values between qubits i and 
j = i + 1 as a function of localization length I. Different 
simbols represents different qubits couples: i — 1 (circles), 
i = 2 (squares), i = 3 (diamonds), i = 4 (triangles), i = 5 
(stars). Arrows point where the exponential decay starts. 
Inset: ^ c as a function of i. The dashed line is an exponential 
fit. 

there is a regime of squared root increasing followed by an 
exponential decay. Notice that here the critical point l c 
where exponential sensitivitivy starts drastically depend 
on the couple of qubits chosen. Indeed, l c ~ 2\ This 
is due to the fact that the number of coefficients in Q 
different from zero, scales as A ^ 2% as the labels of the 
off-diagonal coefficients that are multiplied differ at least 
of 2\ Thus, not negligible off diagonal terms appears 
for greater values of I. In Fig. |3]tbis behavior is shown 
by means of the arrows the predicted points i c (i) = 2 s : 
the value of the localization length corresponding to the 
maximum concurrence saturation value are shown in the 
inset of FigQ] 

We now focus the bipartite entanglement and we char- 
acterizze it by means of the Von Neumann entropy of one 
subsystem. We bipartite the system in two subsystems A 
and B, each composed by tia and tib qubits respectevely, 
we compute 

S A = Trp A \ogp Al (7) 

where pa is the reduced density matrix of the subsys- 
tem A. Notice that, due to the hierarchy introduced 
by the binary coding n = ^ai2 J , the bipartite entan- 
glement displays very different behaviors depending on 
which qubits compose the subsystems. It is necessary 
then to specify both the size and the labels of the qubits 
in each subsystem. We first bipartite the system in qubits 
one to to, varying to (tia — m and ns = n q — to). Then, 
we study the case of subsystem A composed of a single 
qubit (tia — 1 and Ub = n q — l), varying its position. We 
evaluate the Von Neumann entropy after the localization 
time, thus when the wave function of the system can be 
described by Eq.Q and Sa is stationary. 

In FigJSJwe show the saturation level of the entropy of a 



FIG. 5: Von Neumann entropy Sa of a subsystem composed 
by qubits 1 to m for different value of localization length (from 
bottom to top I ~ 2 k , k = 4, 5, 6, 7, 8). 

subsystem of size to for different values of the localization 
length I. We define a critical threshold S c of the entropy 
under which we consider the block of to spins unentan- 
gled with the others. If we choose S c = 1 (straight line 
in Fig[SJ), it is clear that the maximum number of qubits 
entangled increases as log I. It has been shown recently 
[j| that the feasibility of efficient classical simulations of a 
quantum system is conditioned by the presence of bipar- 
tite entanglement and its scaling behavior with respect to 
the number of qubits n. Indeed, if the maximum (with 
respect to the partition) bipartite entanglement in the 
system scales at maximun as log n q , it is possible to per- 
form efficient classical simulations. In our system, we 
may scale the number of qubits keeping fixed the system 
size, thus exploring smaller and smaller scales of the sys- 
tem (T — > 0), or increasing the system size (T = const). 
In the former case I grows exponentially with n q and 
S a ~ n q thus there are no known efficient classical sim- 
ulation methods, while in the latter case I is constant 
and it is thus possible to perform efficient classical sim- 
ulations following 0. This result should not surprise 
as increasing the system size while studying a dynami- 
cally localized system adds almost no information at all, 
as the wave function tail are almost not populated (the 
wave function coefficients decays exponentially). In con- 
clusion, quantum computation of dynamically localized 
systems do improve the efficiency of classical simulations. 
However, the improvement to compute interesting quan- 
tities as the localization length is only quadratic and not 
exponential as shown in [3j. The same arguments apply 
in the simulation of an ergodic quantum chaotic system 
(I = 2 Tlq ). Thus, the simulation of a complex many body 
quantum system, is classically inefficient while it is, at 
least in principle, possible to simulate it efficiently if an 
efficient quantum algorithm exists. 

In Fig. |S|we perform a different analysis, studying the 
Von Neumann entropy of a single qubit with respect to 
the rest of the quantum computer as a function of the 
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FIG. 6: Von Neumann entropy Sa of a subsystem composed 
by a single qubit m for different value of localization length 
(from bottom to top I ~ 2 k , k = 4, 5, 6, 7, 8). 

localization length. The figure shows clearly that the en- 
tanglement depends on the position of the qubit. Again, 
we define a critical threshold under which we consider not 
entangled the qubit. As before, the number of entangled 
qubits scales logarithmically with the localization length. 
Thus, the measure of the reduced density matrix of a sin- 
gle qubit can be used as an efficient method to estimate 
the localization length of a quantum chaotic system. This 
estimate is crucial to investigate complex quantum sys- 
tem efficiently by means of both classical and quantum 



simulations. 

In conclusion we have shown that the concurrence in 
a quantum computer that simulates a dynamically lo- 
calized system is exponentially sensitive to both small 
changes of the Hamiltonian and to the qubits chosen. 
This sensitivity is due to the natural ordering introduced 
on the qubit by the coding of the simulated system. No- 
tice that this is a signature of quantum chaos on a pure 
quantum characteristic with no classical analogue. It 
should be also interesting to compare such sensitivity 
with the cases of classical regular and semi-integrable 
dynamics. Furthermore, the same sensitivity has been 
found recently in the ground state of a single particle 
Anderson model |l7| : these two results reflects the under- 
laying connections between the dynamical and the An- 
derson localization, and a better comprehension of this 
behavior might lay in a more general picture. The results 
on the Von Neumann entropy showed that, the dynami- 
cal regime of the simulated quantum system influence the 
possibility of performing efficient classical simulations. In 
the case of a localized state there are no known meth- 
ods to perform an efficient classical simulation if one is 
exploring smaller scales while increasing the number of 
qubits. A more detailed picture of the dependences of the 
feasibility of efficient quantum simulations depending on 
the dynamical regime of the system will be the objects 
of further studies. 

The author thanks Rosario Fazio and Giuliano Bencnti 
for interesting discussions and a careful reading of the 
manuscript. 
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